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ABSTRACT 

The  Central  Inertial  Guidance  Test  Facility  (CIGTF)  Standard 
Gyrocompasslng  Accuracy  Test  and  associated  data  analysis  procedures  for 
single-degree-of-freedom  gyros  are  presented. 

The  gyrocompasslng  error,  the  error  in  the  north  direction  as 
determined  by  the  gyroscope,  is  obtained  by  orienting  the  gyro  in  four 
specific  orientations  and  measuring  the  output  of  the  gyro  while  in  the 
torque-to-balance  mode.  This  procedure  is  repeated  24  times  for  a  total 
of  25  north  determinations  during  which  the  gyro  is  not  shut  down.  From 
these  data  a  determination  is  made  of  the  average  gyrocompasslng 
error  and  of  the  operating  instability  of  the  gyrocompasslng  error. 

This  "operating"  phase  is  followed  by  a  "shutdown"  phase  consisting 
of  10  north  determinations,  each  separated  by  a  two-hour  gyro  shutdown 
and  two-hour  warmup.  From  these  data  the  shutdown  instability  of  the  gyro 
compassing  error  is  determined.  Although  the  test  and  data  analysis 
procedures  described  in  this  report  are  written  specifically  for  single- 
degree-of-freedom  gyros,  the  theory  can  easily  be  applied  to  two-degree-of 
freedom  gyros. 

This  report  also  contains  a  description  of  the  CIGTF  gyroscope 
laboratory's  two  axis  automated  gyrocompasslng  test  station  and  a  test 
error  analysis. 
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DEFINITIONS  AM)  SYMBOLS 


a  acceleration  component  (2.3) 

a  misalignment  angle  r-Z  the  vertical  fixture  axis  coordinate 

frame  with  respect  r?  the  earth  fixed  coordinate  frame 
measured  about  (Appendix  II) 

b  misalignment  angle  of  the  vertical  fixture  axis  coordinate 

frame  with  respect  to  the  earth  fixed  coordinate  frame  measured 
about  Y-  (Appendix  II) 

c  misalignment  angle  of  the  vertical  fixture  axis  coordinate 

frame  with  respect  to  the  earth  fixed  coordinate  frame  measured 

about  Z_  (Appendix  II) 
b 

D  drift  coefficient  (2.3) 

d  misalignment  angle  of  the  horizontal  fixture  axis  coordinate 

frame  with  respect  to  the  vertical  fixture  axis  coordinate 
frame  measured  about  Xy  (Appendix  II) 

*  degree  (2.3) 

E  (subscript)  earth  fixed  coordinate  frame  (Appendix  II) 

e  misalignment  angle  of  the  horizontal  fixture  axis 

coordinate  frame  with  respect  to  the  vertical  fixture  axis 
coordinate  frame  measured  about  Y^  (Appendix  II) 

F  (subscript)  insensitive  to  acceleration  (2.3) 

f  misalignment  angle  of  the  horizontal  fixture  axis  coordinate 

frame  with  respect  to  the  vertical  fixture  axis  coordinefe 
frame  measured  about  Z^  (Appendix  II) 

G  (subscript)  gyro  axes  coordinate  frame  (Appendix  II) 

g  local  acceleration  of  gravity,  defined  positive  upward  (2.3) 

GSF  gyrocompassing  scale  factor  (4.1) 

H  (subscript)  horizontal  fixture  axis  coordinate  frame  (Appendix  II) 

h  misalignment  angle  of  the  gyro  axes  with  respect  to  the  horizontal 

fixture  axis  coordinate  frame  measured  about  X^  (Appendix  II) 
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Definitions  and  Symbols  (Continued) 
hr  sidereal  hour  (2.3) 

I  (subscript)  input  axis  (Appendix  I) 

I  input  axis  (2.1) 

I  electronically  averaged  value  of  the  torquer  current  over 

a  ten-second  period  (4.1) 

I  average  of  10  ten-second  electronically  averaged  values  of  the 

torquer  current  (2.3) 

i  (subscript)  running  index  for  runs  (5.2) 

i  current  flow  through  the  torquer  (Appendix  I) 

i  misalignment  angle  of  the  gyro  axes  with  respect  to  the 

horizontal  fixture  axis  coordinate  frame  measured  about  Y 
(Appendix  II) 

j  (subscript)  running  index  for  phases  (5.2) 

j  misalignment  angle  of  the  gyro  axes  with  respect  to  the 

horizontal  fixture  axis  coordinate  frame  measured  about 
(Appendix  II) 

k  (subscript)  running  index  for  orientations  (5.1) 

1  (subscript)  running  index  for  the  10  ten-second 

torque-to-balance  current  samples  (5.1) 

ma  milliampere  (2.3) 

N  astronomic  north,  defined  by  the  intersection  of  the 

meridian  plane  with  the  horizontal  plane,  where  the 
meridian  plane  is  defined  by  the  local  gravity  vector 
and  a  line  parallel  to  the  rotation  axis  of  the  earth 
and  the  horizontal  plane  is  defined  as  a  plane 
perpendicular  to  the  local  gravity  vector  (2.2) 

N'  the  north  direction  as  determined  by  the  gyro  (2.2) 

0  (subscript)  output  axis  (2.3) 

0  output  axis  (2.1) 
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Definitions  and  Symbols  (Continued) 


Q  transformation  matrix  (Appendix  II) 

S  (subscript)  spin  axis  (Appendix  I) 

5  spin  axis  (2.1) 

ST  torquer  scale  factor  (2.3) 

s  best  estimate  of  the  run  standard  error  for  a  phase  (5.3) 

V  (subscript)  vertical  fixture  axis  coordinate  frame  (Appendix  II) 

X  axis  of  coordinate  frame  (Appendix  II) 

Y  axis  of  coordinate  frame  (Appendix  II) 

Z  axis  of  coordinate  frame  (Appendix  II) 

B  error  In  180°  rotation  from  position  2  to  position  3 

about  fixture  horizontal  axis  (Appendix  II) 

A  shutdown  instability  of  gyrocompassing  error  (5.3) 

6  operating  instability  of  gyrocompassing  error  (5.3) 

e  error  in  180c  rotation  from  position  1  to  position  2 

about  fixture  vertical  axis  (Appendix  II) 

n  •  error  in  180°  rotation  from  position  3  to  position  4 

about  fixture  vertical  axis  (Appendix  II) 

0  the  gyrocompassing  error:  the  small  angle  measured 

from  astronomic  north  to  the  north  direction  as  determined 
by  the  gyroscope,  positive  as  defined  by  a  rotation 
about  an  axis  directed  vertically  up  (2.2) 

X  local  astronomic  latitude  angle  (Appendix  II) 

£  summation  (5.1) 

o  standard  error  (5.2) 

4>  the  acute  angle  measured  from  a  north-south  line  to  the 

spin  axis,  positive  as  defined  by  a  rotation  about  an 
axis  directed  vertically  up  (2.2) 
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Definitions  and  Symbols  (Continued) 

the  angle  from  north  as  determined  by  the  gyro  to  the 
spin  axis  where  <f> '  =*  $  -  0  (2.2) 

to  rotation  rate  (Appendix  II) 

o)„  earth  rotation  rate  (Appendix  II) 

horizontal  component  of  earth  rotation  rate  (2.3) 
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1. 


INTRODUCTION 


One  of  the  functions  of  a  gyroscope  is  to  provide  an  accurate 
north  reference  for  aircraft  and  missile  guidance  systems.  It  follows 
that,  in  order  to  provide  a  complete  evaluation  of  the  gyroscope,  a 
method  for  determining  the  gyrocompass ing  error,  the  bias  error  in  the 
gyroscope's  determination  of  the  north  direction,  is  necessary.  This 
report  describes  the  method  used  by  the  Central  Inertial  Guidance  Test 
Facility  (CIGTF)  for  evaluating  the  gyrocompass ing  error,  its  operating 
instability,  and  its  shutdown  instability. 

Although  the  test  and  data  analysis  procedures  described  in  this 
report  are  written  specifically  for  single-degree-of-freedom  gyros, 
the  theory  can  easily  be  applied  to  two-degree-of-freedom  gyros.  Also, 
this  method  can  be  used  to  determine  the  effect  of  parameter  variation 
(e.g.,  gyro  rotor  voltage  variation)  and  the  effect  of  gyro  input  axis 
azimuth  variation  on  the  gyrocompassing  accuracy  of  a  gyroscope. 

2.  THEORY 


2.1  Ideal  Case 


Determining  the  gyrocompassing  accuracy  of  a  single-degree- 
of-freedom  gyroscope  involves  orienting  the  gyroscope  so  that  the  input 
axis  (I)  is  constrained  to  the  local  horizontal  plane  and  the  output  axis 
(0)  is  parallel  to  the  local  gravity  vector.  The  gyro  is  connected  in  the 
torque-to-balance  mode,  i.e.,  the  output  of  the  pickoff  is  conditioned 
and  applied  to  the  torquer  with  the  result  that  all  gyro  precessions 
are  "nulled  out".  This  torque-to-balance  current  provides  a  measure  of 
the  total  nulled  precession.  Rotating  the  gyro  input  axis  in  the  horizontal 
plane  (about  the  output  axis-'  until  minimum  "precession"  occurs  would, 
for  the  ideal  single-degree-of-freedom  gyro,  provide  a  means  of 
determining  when  the  input  axis  is  along  an  east-west  line.  When  so 
oriented,  the  input  axis  would  sense  no  component  of  earth  rotation 
rate,  and  as  shown  in  Figure  1,  the  spin  axis  (S)  would  necessarily  be 
in  a  north-south  orientation  and  parallel  to  the  horizontal  component 
of  earth  rate  (to^) . 


UEH  (NORTH) 


FIGURE  1 

MINIMUM  PRECESSION  (Ideal  Case) 


2.2  General  Case 


In  practice,  due  to  misalignment  of  the  gyro  with  respect 
to  the  local  gravity  vector  and  because  of  electrical  and  mechanical 
imperfections  of  the  gyro,  minimum  "precession"  of  the  gyro  will  not 
necessarily  indicate  that  the  spin  axis  is  north.  Even  if  the  gyro 
were  perfectly  aligned  with  respect  to  the  gravity  vector  and  all  internally 
caused  drift  were  compensated  for  based  upon  previous  measurements, 
minimum  "precession"  still  would  not  necessarily  indicate  that  the  spin 
axis  is  precisely  north  because  of  gyro  drift  coefficient  instabilities 
and  coefficient  variations  with  time.  (See  Appendix  I  for  definitions  of 
drift  coefficients.)  The  angle  between  the  spin  axis  and  astronomic 
north*  is  the  angle  $  shown  in  Figure  2. 
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FIGURE  2 

MINIMUM  PRECESSION  (General  Case) 

Physically  $  is  the  sum  of  three  misalignment  angles : 

<!>  =  a  +  d  +  h 

These  three  an:  les  a,  d,  and  h  are  discussed  in  Appendix  IT.  The 
difference  between  true  <p  and  a  computed  <p  (4>')  is  the  gvrocomcassing 
error.  Thus,  it  is  necessary  to  define  an  angle  <f>'  such  that 

<t>'  =  <p  -  6 


*  Astronomic  north  is  defined  by  the  intersection  of  the  meridian  plane 
with  the  horizontal  plane.  The  meridian  plane  is  defined  by  the  local 
gravity  vector  and  a  line  parallel  to  the  rotation  axis  of  th?  earth. 
The  horizontal  plane  is  defined  as  a  plane  nerpendicular  to  .  Jocal 
gravity  vector. 
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where  0  is  the  gyrocorapassing  error  (the  angle  between  north  and  N', 
the  north  direction  as  determined  by  the  gyroscope) . 

To  evaluate  6  in  the  laboratory  the  gyro  is  mounted  on  the 
gyrocompassing  fixture  (section  3)  with  the  spin  axis  approximately 
parallel  to  the  horizontal  axis  of  the  test  fixture  and  initially 
oriented  as  in  Figure  3. 


LABORATORY  ORIENTATION  OF  THE  GYROSCOPE 


By  summing  the  gyro  outputs  from  four  orientations  (section  2.3) 
the  effect  of  an  alignment  error  (angle  h  in  Figure  3)  between 
the  spin  axis  and  the  horizontal  axis  of  the  fixture  is  removed. 
Furthermore,  since  in  the  laboratory  the  horizontal  axis  of  the 
fixture  is  accurately  aligned  to  north  by  optical  methods,  the  angle 
a+d  is  constrained  to  be  zero  degrees.  Therefore  a  solution  for  ' 
in  the  laboratory  is  actually  a  solution  for  the  gyrocompassing 
error  0  . 


2.3  Gyroscope  Orientations 

In  Figure  4  the  four  orientations  used  for  the  test  are 
presented.  These  four  orientations  are  necessary  to  evaluate  0  and 
to  eliminate  the  effects  of  misalignments  (e.g.,  spin  axis  not  parallel 
to  test  fixture  horizontal  axis  and  input  axis  not  horizontal).  A 
discussion  of  these  misalignments  and  their  effects  is  contained  in 
Appendix  II.  If  <p  is  assumed  small,  and  misalignments  (including  h) 
are  neglected,  the  equations  (from  Appendix  1)  relating  the  torque- 
to-halance  current  to  $  in  each  of  the  four  orientations  are 
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FIGURE  4 

GYROCOMPASSING  ORIENTATIONS 
FOR  THE  GYROCOMPASSING  ACCURACY  TEST 
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where 


I1’I2,I3,I4 


“EH 


-  ^eh  +  df  +  Vo 


^UEH  +  DF  +  DOaO 


-  ♦«eh  +  df  +  d0«0 


^EH  +  °F  +  DOaO 


(2) 

(3) 

(4) 


gyro  torquer  scale  factor  ("/sidereal  hour/ma) 

average  torque-to-balance  current  over  10 
ten-second  intervals  In  each  of  the  four 
orientations  (raa) 

horizontal  component  of  earth  rotation 
rate  at  the  Holloman  AFB  latitude 
(12.59562°/sidereal  hr) 


Dp  gyro  drift  which  is  insensitive  to 

acceleration  (°/hr) 

DQa0  gyro  drift  (°/hr)  attributable  to 

^  acceleration  along  the  output  axis 


2. A  Solution  for  0 


If  the  drift  coefficients  in  equations  1  through  A  are  assumed 
to  be  constant  during  the  period  of  the  test,  these  equations  can  be 
combined  to  solve  for  0 . 


'■'hangine  the  si^ns  of  e^uati^ns  1  and  3  and  adding  the  results 
to  equations  2  and  A  gives 


ST  (f2  +  !A 


(L.) 
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Equation  5 
follows  to 


is  independent  of  gyro  drift  and  may  be  rearranged  as 
solve  for  '  (  the  computed  value  of  4>) : 


radians 


(I2  +  h 


(6) 


Since  <j>*  =  <f>  -  6,  and  <j>  is  constrained  to  be  zero  in  the  laboratory, 
-0  can  be  substituted  for  <j> ' .  Then  converting  from  radians  to  arc 
seconds  and  substituting  the  magnitude  of  at  Holloman  AFB  into 
equation  6,  the  gyrocompassing  error  becomes: 


9 sec  '  4093-579  ST  (!1  +  !3  -  l2  -  V 


(7) 


2.5  Gyrocompassing  Scale  Factor 

The  quantity  4093.979  ST  which  has  units  of  arc  seconds/ma 
(or  arc  seconds/mvolt  if  the  torque-to-balance  current  is  measured  by 
inserting  a  precision  resister  in  the  feedback  loop)  is  a  constant 
for  any  one  gyro  and  will  be  called  the  gyrocompassing  scale  factor  (GSF) . 
The  method  used  to  determine  the  GSF  is  described  in  section  4.  Equation 
7  may  now  be  written 


"  GSF  (!1  +  !3  "  h  ~  V 


(8) 


3.  TEST  FIXTURE 


For  the  gyrocompassing  accuracy  test  the  gyro  is  placed  in  a  mount 
which  in  turn  is  mounted  on  the  gyrocompassing  test  fixture.  A  photograph 
of  the  CIGTF  test  station  is  shown  in  Figure  5.  The  gyrocompassing  test 
fixture  is  a  two-axis  test  device  manufactured  by  Goerz  Optical  Company. 

It  is  capable  of  providing  automatic  rotation  about  either  axis  or  about 
both  axes  simultaneously.  The  two  axes,  one  vertical  and  the  other 
horizontal,  are  perpendicular  to  within  five  arc  seconds.  The  mounting 
surface,  which  is  on  the  end  of  the  horizontal  axis,  is  capable  of 
supporting  40  foot-pounds  without  deflecting  the  horizontal  axis  more  than 
0.2  sec  and  the  vertical  axis  more  than  0.5  sec.  The  position  transducer 
on  each  axis  is  an  inductosyn  having  360  stable  nulls  at  one  degree 
intervals  for  360  degrees  of  rotation.  The  fixture  can  thus  be 
positioned  about  either  axis  in  one  degree  increments.  The  positioning 
accuracy  is  better  than  0.5  £ec  with  a  repeatability  of  0.2  sec. 
Positioning  of  the  table  is  accomplished  either  manually  or  by  means  of 
an  automatic  programmer. 
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FIGURE  5 

GYROCOMPASS IMG  TEST  STATION 


4.  TEST  PROCEDURE 


4.1  Determination  of  the  Gyrocompassing  Scale  Factor 

Prior  to  the  start  of  the  test,  the  vertical  test  fixture 
axis  is  aligned  parallel  to  the  local  gravity  vector  by  means 
of  bubble  levels  mounted  on  the  vertical  axis  mounting  surface.  In 
addition,  the  horizontal  axis  is  aligned  parallel  to  the  horizontal 
component  of  earth  rotation  rate  by  optical  means.  Digital  displays 
located  in  the  gyrocompassing  fixture  electronics  racks,  which  indicate 
axis  rotational  position,  should  both  register  zero  degrees.  The  gyro, 
which  is  in  the  torque-to-balance  mode,  is  positioned  with  the  output 
axis  vertical  up  and  the  input  axis  approximately  horizontal  west 
(Orientation  1  in  Figure  4).  After  the  gyro  has  reached  steady  state 
conditions  in  this  orientation,  10  ten-second  electronically  averaged 
values  of  the  gyro  torquer  current,  are  recorded.  The  torque-to-balance 
current  equation  for  this  orientation  is 


Vi  -  -*“eh  +  df  +  Vo  <9) 

The  test  fixture  is  next  rotated  one  degree  in  a  positive 
direction  about  the  vertical  (up)  axis.  After  the  gyro  reaches 
steady  state  conditions  in  this  orientation,  10  ten-second  electronically 
averaged  values  of  the  gyro  torquer  current,  I',  are  recorded.  In  this 
orientation  the  torque-to-balance  current  equation  is 


Vi  -  -  ('♦  +  .01745)  -n,  +  Bp  +  v0 


(10) 


where 


1°  =  .01745  radians 


and  1^ 


is  the  average  torque-to-balance  current  of  the  ten  I|. 


Then 


V1!  - 


)  ■  .01745  to 


EH 


(ID 


8 


When  both  sides  of  the  equation  are  divided  by  4to„  (I.. 
the  result  is 


9 


-♦-0.0136., . .  radians/ma 

<!!  -  !i> 


(12) 


Converting  radians  to  sec,  and  substituting  the  value  for  the 
horizontal  component  of  earth  rotation  rate,  yields 


4093.979  ST  - 


900 


sec/ma 


(13) 


Thus  the  gyrocompassing  scale  factor  GSF  is  given  by 


GSF 


900 


sec/ma 


(14) 


4.2  Four  Orientation  Test  Sequence 

The  gyrocompassing  test  sequence  which  has  been  programmed 
on  paper  tape  is  initiated  by  the  test  fixture  automatic  programmer. 
The  sequence  of  operations  that  are  produced  by  the  programmer  are  as 
follows : 


4.2.1  With  the  spin  axis  parallel  to  the  horizontal  fixture 
axis,  the  gyro  is  rotated  about  the  fixture  horizontal  and/or  vertical 
axes  until  the  gyro  is  in  Orientation  1  (0  vertical  up,  I  horizontal 
west).  This  orientation  is  maintained  until  the  gyro  has  settled  and 
10  ten-second  electronically  averaged  values  of  the  gyro  torquer  current 
are  recorded.  The  values  are  automatically  printed  and  punched  on  paper 
tape. 


4.2.2  The  gyro  is  rotated  180°  about  the  fixture  vertical 
axis  to  place  the  gyro  in  Orientation  2  (0  vertical  up,  I  horizontal  east). 
This  orientation  is  maintained  until  the  gyro  has  settled  and  10  ten-second 
electronically  averaged  values  of  the  gyro  torquer  current  are  recordec. 


% 
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4.2.3  Next  the  gyro  is  rotated  180°  about  the  fixture 
horizontal  axis  to  place  the  gyro  in  Orientation  3  (0  vertical  down,  I 
horizontal  west).  This  orientation  is  maintained  until  the  gyro  has 
settled  and  10  ten-second  electronically  averaged  values  of  the  gyro 
torquer  current  are  recorded. 

4.2.4  The  gyro  is  rotated  180°  about  the  fixture  vertical 
axis  to  place  the  gyro  in  Orientation  4  (0  vertical  down,  I  horizontal 
east).  This  orientation  is  maintained  until  the  gyro  has  settled 

and  10  ten-second  electronically  averaged  values  of  the  gyro  torquer 
current  are  recorded. 

4.3  Operating  Phase 

The  sequence  of  commands  and  the  resulting  data  obtained  from 
the  gyro  in  the  four  orientations  described  in  section  4.2  are  necessary 
for  one  north  determination  or  one  "run".  This  procedure  is  repeated  24 
times  for  a  total  of  25  runs  during  which  the  gyro  is  not  shut  down.  This 
comprises  the  "operating  phase"  of  the  gyrocompassing  accuracy  test. 

4.4  Shutdown  Phase 


The  operating  phase  is  followed  by  a  gyro  shutdown  of 
approximately  two  hours.  After  the  two  hour  shutdown  the  gyro  is 
warmed  up  for  two  hours  and  one  north  determination  is  made.  This 
procedure  of  shutdown,  warmup,  and  north  determination  is  repeated 
nine  times.  These  ten  north  determinations,  each  separated  by  a 
gyro  shutdown  and  warmup  period,  comprise  the  "shutdown  phase" 
of  the  gyrocompassing  accuracy  test. 

5.  DATA  ANALYSIS 


5 . 1  Orientation  Values 


5.1.1  The  mean  value  of  the  10  ten-second  torque-to-balance 
current  samples  in  orientation  k  is : 


I 


k 


10 


I 

£-1 


(15) 


where 

k  =  1,  2,  3,  and  4 


10 


} 
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5.1.2 

orientation  k  is: 


The  variance  of  the  10  current  samples  in 


10 


l 

Jl=l 


(16) 


5.1.3  The  variance  of  the  mean  of  the  10  current  samples 
in  orientation  k  is : 


(17) 


5.2  Run  Values 

5.2.1  The  gyrocompassing  error  in  arc  seconds  for  run 
i  of  phase  j  is: 


6 


ij 


I5SF  (I1  +  f3  -  i2  -  I4> 


(18) 


where 

jsl  for  the  operating  phase  and  j=2  for  the  shutdown  phase. 

5.2.2  The  best  estimate  of  the  standard  error  of  0  for 
run  i  of  phase  j  is : 


r 

f, 

i 


5 . 3  Phase  Values 


5.3.1  The  gyrocompassing  error  for  phase  i  is: 


6 


j 


n 


I 

i=l 


n 


(20) 


where  n  is  the  number  of  runs  in  the  phase. 

5.3.2  The  best  estimate  of  the  standard  error  ui  0.  is: 


o 


1 

n 


(21) 


except  for  the  case  described  in  paragraph  5.4.2. 


5.3.3  The  best  overall  estimate  (RMS)  of  the  run 
standard  error  for  the  n  runs  constituting  phase  j  is: 


which  by  equation  (21)  becomes 


(22) 


s 


j 


(23) 
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5.3.4  The  operating  instability  is  the  best  estimate 
of  the  standard  deviation  of  the  0^  from  the  operating  phase,  i.e.. 


6 


(24) 


where  j«l. 


5.3.5  The  shutdown  instability  is  the  best  estimate  of 
the  standard  deviation  of  the  0  from  the  shutdown  phase,  i.e., 


(25) 


where  j=2. 

5.4  Significance  Tests 

5.4.1  A  Student  "t"  test  is  performed  to  determine  if  0 
is  significantly  different  from  zero  at  the  5*  significance  level.  For 
a  large  number  of  degrees  of  freedom 


6i  I  ;>  1.96 

implies  that  the  average  gyrocompassing  error  for  phase  j  is  significant. 
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5.4.2  A  Snedecor  F  Test  is  performed  to  determine  if 
the  operating  instability  is  significantly  greater  than  the  RMS  of  the 
run  standard  errors  at  the  5%  significance  level.  With  24  degrees 
of  freedom  for  the  numerator  and  900  degrees  of  freedom  for  the 
denominator, 


A2 

*  1.53 


implies  that  the  variation  introduced  by  repositioning  the  gyro  is 
significantly  greater  than  that  attributable  to  noise,  and  the  best 
estimate  of  the  standard  error  of  0,  is  then  —  rather  than  on  . 

J  ^  ej 

In  this  case  the  "t"  test  described  in  5.4.1  is  corrected  accordingly. 

5.4.3  A  Snedecor  F  Test  is  performed  to  determine  if  the 
shutdown  instability  is  significantly  greater  than  the  F>tS  of  the 
run  standard  errors  of  the  shutdown  phase  at  the  5%  significance  level. 
With  9  degrees  of  freedom  for  the  numerator  and  360  degrees  of  freedom 
for  the  denominator,  then 


> 


1.91 


implies  that  the  variation  introduced  by  shutdown  and  by  repositioning 
the  gyro  is  significantly  greater  than  that  attributable  to  noise. 

5.4.4  A  Snedecor  F  Test  is  performed  to  determine  if 
the  shutdown  instability  is  significantly  greater  than  the  operating 
instability  at  the  57.  significance  level.  With  9  degrees  of  freedom 
for  the  numerator  and  24  degrees  of  freedom  for  the  denominator,  then 


-  2.30 


implies  that  the  variation  introduced  by  gyro  shutdown  between  runs 
is  significantly  greater  than  that  attributable  to  repositioning  the 
gyro. 
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APPENDIX  I 

S INGLE-DEGREE-OF-FREEDOM  GYRO  PERFORMANCE  MODEL 

The  assumed  performance  model  (Reference  1:2-4)  for  a  single- 
degree-of-freedom  gyro  is 


STi  “  °F  +  Diai  +  Doao  +  Dsas  +  Dnai  +  Dssas  +  Dioaiao 


+  Disaias  +  Dosaoas  +  WI 


where 

i  current  flow  through  the  torquer  (ma) 

ST  torquer  scale  factor  (°/hr/ma) 

a^.  acceleration,  with  respect  to  inertial  space, 

1  of  the  gyro  case  along  the  input  axis  ig) 

aQ  acceleration,  with  respect  to  inertial  space, 

of  the  gyro  case  along  the  output  axis  (g) 

ag  acceleration,  with  respect  to  inertial  space, 

of  the  gyro  case  along  the  spin  axis  (g) . 

Dp  gyro  drift  (°/hr)  which  is  insensitive  to 

acceleration 

D^a^  gyre  drift  (°/hr)  attributable  to  acceleration 

along  the  input  axis,  where  (°/hr/g)  is  a 
drift  coefficient 

Dq3q  gyro  drift  (°/h'r)  attributable  to  acceleration 

along  the  output  axis,  where  (°/hr/g)  is  a  drift 
coefficient 

D„a  gyro  drift  (°/hr)  attributable  to  acceleration  along 

the  spin  axis,  where  D^  (°/hr/g)  is  a  drift  coefficient 

2 

D^a^  gyro  drift  (°/hr)  attributable  to  the  square  of  2 

1  '  acceleration  along  the  input  axis,  where  (°/hr/g  ) 

is  a  drift  coefficient 


* 


Dssas 


Dioaiao 


Disaias 


DosVs 


WI 

hr 


gyro  drift  (°/hr)  tributable  to  the  square  of 
acceleration  along  .ne  spin  axis,  where  D  (°/hr/g^) 
is  a  drift  coefficient 


gyro  drift  (°/hr)  attributable  to  the  product  of 
accelerations  along  the  input  axis  and  output  axis 
where  (°/hr/g^j  is  a  drift  coefficient 


gyro  drift  (°/hr)  attributable  to  the  product  of 
accelerations  along  the  input  axis  and  spin  axis, 
where  D^s  (°/hr/g2)  is  a  drift  coefficient 

gyro  drift  (°/hr)  attributable  to  the  product 
of  accelerations  along  the  output  axis  and  spin 
axis,  where  D0g(6/hr/g2)  is  a  drift  coefficient 


angular  velocity,  with  respect  to  inertial  space, 
of  the  gyro  case  about  the  input  axis  (°/hr) 


sidereal  hour 


g 


local  acceleration  of  gravity,  defined  positive 
upward 
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APPENDIX  II 


EFFECT  OF  MISALIGNMENT  ANGLES 
IN  THE  GYROCOMPASSING  ACCURACY  TEST 

1.  INTRODUCTION 

In  the  laboratory  gyrocompassing  accuracy  test  the  effect  of 
misalignment  angles  must  be  considered  and,  where  possible,  eliminated, 
fhe  identification  and  elimination  of  the  misalignment  angles  are 
discussed  here. 

2.  DERIVATION  OF  DRIFT  RATE  DUE  TO  MISALIGNMENTS 
2.1  Transformation  of  Coordinates 

In  order  to  derive  an  expression  for  all  possible  misalignment 
angles,  three  transformations  of  coordinates  for  each  of  the  four 
test  orientations  are  necessary.  The  transformations  will  be  accomplished 
in  the  following  sequence:  the  earth  fixed  coordinates  (defined  by  the 
local  gravity  vector  and  astronomic  north)  will  be  transformed  into  the 
vertical  fixture  axis  coordinate  frame;  the  vertical  fixture  axis 
coordinate  frame  will  be  transformed  into  the  horizontal  fixture  axis 
coordinate  frame;  and  finally  the  horizontal  fixture  axis  coordinate 
frame  will  be  transformed  into  the  gyro  input,  output,  and  spin  axes 
coordinate  frame. 


Y 

Z_  (WEST) 

h. 

FIGURE  II-l 

EARTH  FIXED  COORDINATE  FRAME 


% 


In  the  earth  fixed  coordinate  frame  shown  above,  the  local  gravity 

vector  g  is  parallel  to  the  X„  axis,  Y_  is  astronomic  north,  and  Z 
,  fc,  fc.  i- 


Angles  a,  b,  and  c  are  associated  with  the  earth  fixed 
to  vertical  fixture  axis  transformation  of  coordinates  in  the 
following  manner: 


cos  b  0 

0  1 

sin  b  0 


— 

“  — 

-sin  b 

XE 

0 

ye 

cos  b 

ZE 

_  _ 

1 1 —2 


t 


where  is  the  vertical  fixture  axis  and  Yy  and  Z  are  in  a  plane 
perpendicular  to  it.  Using  small  angle  approximations 


sin  x  =  x 

r 

i  cos  x  =  1 

and  neglecting  second  degree  terms,  the  vertical  fixture  axis  coordinate 
:  frame  may  be  written  in  terms  of  the  earth  fixed  coordinate  frame: 

* 


> 

*v 

1  c  0 

1  0  -b 

10  0 

XE 

Y„ 

a 

-c  1  0 

0  10 

Ola 

Y_ 

% 

r 

V 

E 

zv 

0  0  1 

b  0  1 

0  -a  1 

ZE 

_  m 

_  m 

m  _ 

„  _ 

_  _ 

or 
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*v 

1  c  -b 

XE 

Ytr 

S3 

-c  1  a 

Y_, 

V 

E 

z„ 

b  -a  1 

Zr. 

V 

E 

- 

Therefore  the  transformation  matrix  is  defined: 


1  c  -b 

— c  1  a 

b  -a  1 


Similarly,  angles  d,  e,  and  f  are  associated  with  the  transformation 
of  the  gyrocompassing  fixture  vertical  axis  coordinates  to  the 
gyrocompa8sing  fixture  horizontal  axis  coordinates.  YH  is  the 
horizontal  fixture  axis  with  and  perpendicular  to  it.  The 
corresponding  transformation  matrix  is  where 


i 

i 


1 


f  -e 


i 


-f 


1 


d 


e 


-d 


1 


Angles  h,  i,  and  j  are  associated  with  the  transformation  of  the 
gyrocompassing  fixture  horizontal  axis  coordinates  to  the  gvro 
coordinate  frame.  The  transformation  matrix  is  Q„u,  where 


II-4 


/ 


Therefore  the  transformation  matrix  for  earth  fixed  to  gyro  coordinate 
frame  is 


1  (c+f+j)  -(b+e+i) 


QGE  "  QGH  qhv  qve 


-(c+f+j)  1 

(b+e+i)  -(a+d+h) 


(a+d+h) 

1 

. 


2.2  Orientation  One 


In  orientation  one  of  the  gyrocompassing  accuracy  test  the 
output  axis  (Xg)  is  up,  the  input  axis  (Z^)  is  west,  and  the  spin 
axis  (Y^)  is  notch.  The  resulting  components  of  gravity  along  the 
three  axes  in  units  of  g  are 


1 

CD 

o 

"  1  ' 

as 

0 

SI 

0 

_  - 

_  _ 

where  Qr,_  is  the  Q„_.  defined  in  section  2.1  of  this  appendix. 
Therefore 

ao  15  1 

ag  =>  -  (c+f+j) 

*=  +  (b+e+i) 
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The  components  of  earth  rate  along  these  three  axes  are 


r  n 

— 

o 

3 

to£  sin 

A 

“s 

Q°E1 

i0_  cos 
E 

A 

ui 

0 

- 

where  is  earth  rotation  rate  and  A  is  the  astronomic  latitude  angle. 
The  only  component  sensed  by  the  gyro  is  : 


a)T  =  (b+e+i)  a)„  sin  A  -  (a+d+h)  cos  A 

It,  L 


When  these  quantities  are  substituted  into  the  assumed  performance 
model,  (Appendix  I),  the  resulting  drift  equation  for  orientation  1  is: 


^1*1  “  Df  +  ^+e+i)  Dx  +  °0  "  +  (k+e+i)  Djq 

-  (c+f+j)  D  +  (b+e+i)  u>_,  sin  A  -  (a+d+h)  w  cos  A 

Uo  t  t 


2.3  Orientation  Two 


To  arrive  at  orientation  2,  a  rotation  of  180°  +  e 
about  the  gyrocompassing  fixture  vertical  axis  is  required,  where  e 
is  a  small  error  in  the  180°  rotation.  Therefore 


QGHQHV 


10  0 
0  -1  -e 


0 


€ 


-1 


Thus 


t 


1  (c-f-j)  (-b+e+i) 

Qge  *  (c-f-j)  -1  -(a+d+h+e) 

(-b+e+i)  (a+d+h+e)  -1 

and  the  components  of  acceleration  along  the  three  gyro  axes  are: 

ao  =  1 

ag  “  c-f-j 

3I  =  -b+e+i 

The  earth  rate  sensed  by  the  gyro  is 

wT  3  (-b+e+i)  sin  X  +  (a+d+h+e)  u>r  cos  X 

"Tie  resulting  drift  equation  for  orientation  2  is: 

St!2  -  Dp  +  (-b+e+i)  Dr  +  DQ  +  (c-f-j)  Dg  +  (-b+e+i)  DI0 

+  (c-f-j)  D  +  (-b+e+i)  u)  sin  X  +  (a+d+h+e)w  cos  X 

UO  L  L 
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sf 
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Orientation  Three 


2.4 


To  arrive  at  orientation  3  it  is  necessary  to  rotate  180°  +  8 
about  the  gyrocompassing  fixture  horizontal  axis  Y^,  where  8  is 
an  error  in  the  180"'  rotation.  Therefore 


-1 

0 

8 

1 

0 

~1 

0 

■  qgh 

0 

1 

0 

qhv 

0 

-1 

-e 

-8 

0 

-1 

0 

e 

-1 

_ 

. 

-1 

(-c+f-j) 

(b-e-i-8) 

a 

(c-f+j) 

-1 

(-a-d+h-e) 

(b-e-i-8) 

(-a-d+h  -t ) 

1 

The  components  of  acceleration  along 

the  gyvo 

are 

as  =  c-f+j 
=  b-e-i-8 


an<i  the  earth  rate  sensed  by  the  gyro  is 

wT  «  (b-e-i-8)  w_  sin  X  -  (a+d-h+e)  u)„  cos  X 

I  £  £ 


The  resulting  drift  equation  for  orientation  3  is: 


STI3  =  Dp  +  (b-e-i-8)  Dj  -  DQ  +  (c-f+j)  Dg  -  (b-e-i-8)  DI0 


(c-f+j)  D„„  +  (b-e-i-S)w,.  sin  X  -  (a+d-h+e)  uj„  cos  X 
Ub  £  fc. 
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Orientation  Fo 
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To  arrive  at  orientation  four,  it  is  necessary  to  rotate  180°+n 
about  the  gyrocompassing  fixture  vertical  axis,  Then 


r*  -i 

-10  8 

10  0 

10  0 

0  10 

qhv 

0  -l  -n 

0  -1  -e 

-8  0  -1 

0  n  -1 

0  £  -1 

—  «. 

-1 


(-c-f+j) 


(-c-f+j) 


(b+e+i+8) 


(a+d-h+e+n ) 


jj[-b-e-i-6)  (a+d-h+tJ-n)  -1 


The  components  of  acceleration  along  the  gyro  axes  are 

a0  =  *  1 
as  =  -(c+f-j) 

=  -(b+e+i+B) 

and  the  earth  rate  sensed  by  the  gyro  is 

=  -  (b+e+i+8)  u)g  sin  A  +  (a+d-h+e+n)  u>g  cos  A 

Therefore  the  drift  equation  for  orientation  four  is: 

St!4  -  Dp  -  (b+e+i+8)  -  DQ  -  (c+f-j)  Dg  +  (b+e+i+6)  DJ0 

+  (c+f-j)  D  -  (b+e-t-i+8)  w_  sin  A -4  (a+d-h+c+n)  cos  A 

vO  it  E* 


t 


2 . 6  Resulting  Drift  Equation 

When  the  drift  equations  for  orientations  1  and  3  are 
subtracted  from  the  sun  of  those  for  orientations  2  and  4,  the  following 
equation  results. 


ST(I2  +  h 


X1  “ 


V 


-  4bDI  + 


4c  D 


OS 


-  4b^_  sin 
t 


+  4  ( a+d ) 


cos 


+  (3c  +  n)  ,ji  cos  ' 
E 


(26) 


The  angle  (a+d)  is  the  angle  between  north  and  the  horizontal  fixture 
axis,  measured  positive  about  an  axis  directed  vertically  un.  Tc 
solve  for  its  computed  value  (a+d)'  equation  26  is  rearranged  as  follows 


(a+d) '  a 


ST(I2  +  I4  -  !x  -  !j)  b(*E  sin  •  +  Dt) 

4  ‘jj  cos  1  cos  - 

E  E 


cD0?  ( 3> +r ) 
.^cos'  4 


(27) 


The  angles  b  and  c  are  the  realignment  angle:  between 
the  vertical  axis  of  the  gyrocompassing  test  fixture  and  the  local 
gravity  vector  in  the  east-west  plane  and  the  north-south  piane^ 
respectively.  These  misalignments  are  minimized  (less  than  1  set) 
by  the  use  of  bubble  levels  prior  to  gyrocompassing.  The  angles  .  and  - 
are  errors  in  the  180°  rotation  between  orientations.  Ihe  positioning 
accuracy  of  the  gyrocompassing  test  fixture  Is  0.5  £ec  or  better. 

Since  these  misalignment  effects  are  verv  small,  thev  are  neglected, 
thus 


(a+d)  ' 


VX2  +  X4 


it 


COS  ‘ 


(28) 


From  section  2.2 


*  =  fi  -  6 

and 

f  =  a+d+h 


(  24  ) 


(  30  i 
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Since  h  is  not  computed  equation  30  may  be  written 


$ '  «  ( a+d ) 1  +  h 


(31) 


where  <j>'  is  the  computed  value  of  $  and  (a+d)'  is  the  com-uteo  value 
of  (a+d).  Substituting  equations  30  and  31  into  equator.  29  an'’ 
subtracting  h  from  both  sides  of  the  equation  yields 


(a*d)’  =  (a+d)  -  9 


In  the  laboratory  test  the  angle  a+d  is  constrained  to  be  zero. 

Thus  a  solution  for  (a+d)'  becomes  a  solution  for  -9.  From  equation 
28  the  gyrocompass ing  error  then  becomes 


radians 


h  »i  *  ;3  -  h  -  v 

4  w_  cos  * 


'.'Then  radians  are  converted  to  arc  seconds  and  the  magnitude  of  cos 
at  Holloman  AFB  is  substituted  into  this  equation,  the  result  is 
equation  7  of  section  2.4. 

Therefore,  the  effects  of  mi  sail gnTents  ’H.ch  can  not  be 
readilv  estimated  or  minimized,  such  as  the  mi  sal  i  m<.  nt  of  the 
spin  axis  with  respect  to  the  horizontal  fixture  axi1:  (angle  h)  and 
the  misalignment  of  the  gvro  output  axis  with  respect  to  the 
vertical  fixture  axis  (angles  i,  1,  e,  and  f),  are  removed  from  the 
determination  of  the  gyrocompasslnv  erroi  bv  the  f our- or i entat i on  test. 
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»3  ABSTRACT 


The  Central  Inertial  Guidance  Test  Facility  (CIGTF)  Standard  Gyrocompassing 
Accuracy  Test  and  issociated  data  analysis  procedures  for  single-degree-of-freedom 
gyros  are  presented. 

The  gyrocompassing  error,  the  error  in  the  north  direction  as  determined  by  the 
gyroscope,  is  obtained  by  orienting  the  gyro  in  four  specific  orientations  and 
measuring  the  output  of  the  gyro  while  in  the  torque-to-balance  mode.  This  procedure 
is  repeated  24  times  for  a  total  of  25  north  determinations  during  which  the  gyro  is 
not  shut  down.  From  these  data  a  determination  is  made  of  the  average  gyrocompassing 
error  and  of  the  operating  instability  of  the  gyrocompassing  error.  This  "operating" 
phase  is  followed  by  a  "shutdown"  phase  consisting  of  10  north  determinations,  each 
separated  bv  a  two-hour  gyro  shutdown  and  two-hour  warmup.  From  these  data  the 
shutdown  instability  of  the  gyrocompassing  error  is  determined.  Although  the  test 
and  data  analysis  procedures  described  in  this  report  are  written  specifically  for 
sing’.e-degree-of-f reedom  gyros,  the  theory  can  easily  be  applied  to  two-degree-of- 
freedom  gyros. 

This  report  also  contains  a  description  of  the  CIGTF  gyroscope  laboratory's 
two  axis  automated  gyrocompassing  test  station  and  a  test  error  analysis. 
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